Two-channel Kondo physics in odd impurity chains 
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We study odd-membered chains of spin-i impurities, with each end connected to its own metallic 
lead. For antiferromagnetic exchange coupling, universal two-channel Kondo (2CK) physics is shown 
to arise at low energies. Two overscreening mechanisms are found to occur depending on coupling 
strength, with distinct signatures in physical properties. For strong inter-impurity coupling, a 
residual chain spin-* moment experiences a renormalized effective coupling to the leads; while in 
the weak-coupling regime, Kondo coupling is mediated via incipient single-channel Kondo singlet 
formation. We also investigate models where the leads are tunnel-coupled to the impurity chain, 
permitting variable dot filling under applied gate voltages. Effective low-energy models for each 
regime of filling are derived, and for even-fillings where the chain ground state is a spin singlet, an 
orbital 2CK effect is found to be operative. Provided mirror symmetry is preserved, 2CK physics is 
shown to be wholly robust to variable dot filling; in particular the single-particle spectrum at the 
Fermi level, and hence the low-temperature zero-bias conductance, is always pinned to half-unitarity. 
We derive a Friedel-Luttinger sum rule and from it show that, in contrast to a Fermi liquid, the 
Luttinger integral is non-zero and determined solely by the 'excess' dot charge as controlled by gate 
voltage. The relevance of the work to real quantum dot devices, where inter-lead charge-transfer 
processes fatal to 2CK physics are present, is also discussed. Physical arguments and numerical 
renormalization group techniques are used to obtain a detailed understanding of these problems. 

PACS numbers: 71.27.+a, 72.15.Qm, 73.63.Kv, 71.10.Hf 



I. INTRODUCTION 

Non-Fermi liquid (NFL) behavior is famously realised 
in the two-channel Kondo (2CK) model, 1 in which a sin- 
gle spin- ^ impurity is exchange-coupled to two equivalent 
but independent metallic conduction bands. Its fascina- 
tion for theorists is evident in the wide range of tech- 
niques applied to the model, including notably the Bcthc 
ansatz, 2-6 numerical renormalization group 6-10 and con- 
formal field theory 10-13 (for a review, see Ref. 14). Such 
methods have elucidated key aspects of the NFL state 
arising from 'overscreening' of the impurity spin below 
the low-energy 2CK scale, 1 T%P K ; including exotic phys- 
ical properties such as a residual entropy of \ hi(2), and a 
logarithmically divergent low-temperature magnetic sus- 
ceptibility, which are symptomatic of the frustration in- 
herent when two conduction channels compete to screen 
the impurity local moment. 14 NFL scaling of the conduc- 
tance has also been predicted theoretically 13-17 , and mea- 
sured experimentally 18 in quantum dots, with its square- 
root temperature dependence being a characteristic sig- 
nature of the 2CK phase. 

This is all of course in marked contrast to standard 
Fermi liquid (FL) behavior, 19 arising for example in the 
single-channel spin-i Kondo or Anderson models 19 (re- 
alised in practice in e.g. ultrasmall quantum dots 20 ' 21 ). 
Here, the dot spin is completely screened by a single 
bath of conduction electrons. The impurity entropy is in 
consequence quenched on the lowest energy scales, and 
the susceptibility is constant. 19 Such systems are char- 
acterised by a unitarity zero-bias conductance, with a 



quadratic temperature dependence at low- energies. 2 2 

But the NFL physics of the 2CK model is delicate: 
finite channel-asymmetry and/or inter-channel charge 
transfer ultimately drive any real system out of the NFL 
regime, 11,14 to a FL ground state. The exquisite tun- 
ability of quantum dot devices 23 allows manipulation of 
such perturbations; indeed couplings can be fine-tuned 
via application of gate voltages to effectively eliminate 
channel-asymmetry. However, tunnel-coupling in such 
systems must result in some degree of charge transfer 
between the metallic 'leads'. This is of course responsi- 
ble for the predominance of single-channel Kondo physics 
in real quantum dot systems. 

Suppressing inter-channel charge transfer allows for 
the emergence of 2CK physics at intermediate tempera- 
tures/energies, although the instability of the 2CK fixed 
point to charge-transfer means that an incipient NFL 
state forming at T ~ Tjf K is subsequently destroyed 
below a FL crossover scale Tfl • Observation of NFL be- 
havior at higher temperatures thus depends on a clear 
separation of scales, Tfl "C Tff K . This was achieved 
recently 18 through use of an interacting lead tuned to 
the Coulomb blockade regime, and to date it is the only 
unambiguous experimental demonstration of the 2CK ef- 
fect. Alternatively, sequential tunneling through sev- 
eral coupled dots should suppress charge transfer be- 
tween leads, 24 and this could be exploited to access 2CK 
physics; although the interplay between spin and or- 
bital degrees of freedom in coupled quantum dot sys- 
tems can also generate new physics, as well known both 
theoretically 24-46 and experimentally. 18 ' 47-56 
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In light of the above, we here consider odd-mcmbcred 
coupled quantum dot chains (each end of which is con- 
nected to its own metallic lead), and demonstrate that 
2CK physics is indeed generally accessible in these sys- 
tems. In models where the couplings are of pure exchange 
type, we show that the low-energy behavior is described 
by the channel-asymmetric two-channel Kondo model, 14 
with pristine 2CK physics surviving down to the lowest 
energy scales in the mirror-symmetric systems of most 
interest. Such models are considered in Sec. II, where 
analytic predictions are confirmed and supplemented by 
use of Wilson's numerical renormalization group (NRG) 
technique 57-61 (for a recent review see Ref. 59). In par- 
ticular, universal scaling of thermodynamic and dynamic 
properties is demonstrated for odd chains of different 
length, and for systems with different coupling strengths. 
While the underlying 2CK physics of such systems is 
shown to be robust for finite antiferromagnetic exchange 
coupling, we find that the mechanism of overscreening 
differs according to whether the inter-impurity coupling 
is strong or weak. In the former case (Sec. II A), it is 
the single lowest spin-^ state of the chain which effec- 
tively couples to and is overscreened by the leads; while 
for weak inter-impurity coupling (Sec. II B), two-channel 
Kondo coupling is mediated via incipient sm^Ze-channel 
Kondo singlets. Clear signatures of the latter are evident 
in the behavior of the frequency-dependent t-matrix, re- 
sults for which are presented and analysed. 

In order to investigate conductance across two-channel 
coupled quantum dot chains, we study in Sec. Ill a re- 
lated class of models in which the terminal impurities are 
Anderson-like quantum levels, tunnel-coupled to their re- 
spective metallic leads (although with inter-lead charge 
transfer still precluded by inter-impurity exchange cou- 
plings). In the mirror-symmetric systems considered, we 
derive effective low-energy 2CK models valid for each 
regime of electron- filling. Even-occupation filling regimes 
- where the chain ground state is a spin-singlet - are 
found in particular to exhibit an orbital 2CK effect, with 
spin playing the role of a channel index. In consequence, 
2CK physics is found to be robust throughout all regimes 
of electron-filling induced by changes in gate voltage. 

Single-particle dynamics for such systems are then con- 
sidered, and hence conductance (Sec. Ill CI), the S- 
matrix and associated phase shifts (Sec. HID); again 
highlighting the universality arising at low-energies. It 
is found in particular that the Fermi level value of the 
T = single-particle spectrum - and hence the zero-bias 
conductance - is pinned to a half-unitary value, irrespec- 
tive of electron-filling. A Fricdcl-Luttinger sum rule 62 is 
then derived in Sec. HIE, relating the Fermi level value 
of the spectrum to the 'excess' charge due to the quan- 
tum dot chain, and the Luttinger integral. 63,64 By virtue 
of the spectral pinning, the sum rule relates directly the 
Luttinger integral to the excess charge/dot filling; in con- 
trast to a Fermi liquid, where it is the Luttinger integral 
that is ubiquitously 'pinned' (to zero), 63 ' 64 and the dot 
filling then determines the value of the spectrum at the 



Fermi level. 19,65 

Finally, in Sec. IV we consider briefly the applicability 
of our findings to real coupled quantum dot devices. We 
argue that the effective low-energy model describing such 
systems is generically a 2CK model with both channel- 
asymmetry and inter-lead cotunneling charge transfer. 
Via a suitable basis transformation, one obtains a model 
in which charge transfer between even and odd chan- 
nels is eliminated, whence the underlying behavior is 
readily understood in terms of that of a pure channel- 
asymmetric 2CK model. The spectrum/t-matrix in a 
given physical channel is however related via this trans- 
formation to a combination of t-matrices in even and odd 
channels. In the mirror-symmetric case sought experi- 
mentally, this leads to the striking conclusion that for 
sufficiently small but non-vanishing cotunneling charge 
transfer, the crossover out of the NFL regime is not in 
fact apparent in conductance measurements, despite the 
ultimate low-energy physics being that of a Fermi liquid. 

II. 2CK HEISENBERG CHAINS 

We consider a chain of N c coupled spin- 1 impurities, 
each end of which is also coupled to its own metal- 
lic lead, as illustrated in Fig. 1. To investigate 2CK 
physics, we study explicitly in this section a system 
of exchange-coupled spin-| impurities, where inter-lead 
charge transfer is eliminated from the outset. For an 
impurity chain of length N c , the full Hamiltonian is thus 
H»-=H L + H?-+H»*. Hereff i = Ea,k,^ k at kff a ak(T 
refers to the two equivalent non-interacting metallic leads 
(a = L/R), and 

N c -1 

H c c =jJ2Si- Sm (la) 

i=l 

H^ c —Jkl^i ■ slo + Jkr&Nc ■ s_ro (lb) 

where Sj is a spin- i operator for impurity i, and s Q o is 
the spin density of lead a = L(R) at impurity 1(N C ): 

<j,(j' 

k 

with <r the Pauli matrices and / a0cr the creation operator 
for the 'O'-orbital of the a = L/R Wilson chain (N — > oo 
is the number of orbitals/k-states in the chain). 

In the following, we focus on odd-7V c chains, with an- 
tiferromagnetic (AF) exchange couplings for both the 
intra-chain Heisenberg exchange, J > 0, and the Kondo 
exchanges, Jku > 0. [We do not consider here even- 
N c chains, where the generic low-energy physics is quite 
different and in essence that of the two-impurity Kondo 
model]. We also consider both the channel symmetric 
case, Jkr = Jkl, as well as the more general case where 
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FIG. 1: Schematic illustration of odd-membered impurity 
chains, each end of which is coupled to its own metallic lead. 

channel asymmetry is present, Jkr. ^ Jkl- The sim- 
plest member of the family, N c = 1, is of course the 
classic single-spin 2CK model, 1 ' 14 while variants of the 
N c = 3 trimer have also been considered previously. 38-45 
We show below that an effective single-spin 2CK model 
results in all cases, whence universal 2CK physics is 
expected on the lowest energy scales in the channel- 
symmetric case; albeit that the mechanism by which the 
effective two-channel coupling arises is rather different in 
the strong and weak inter-impurity coupling regimes. 

A. Strong inter-impurity coupling 



In the absence of a magnetic field, t / 4- spin sym- 
metry implies Pf±\N c ; S z ) = i\N c ; —S z ), where P-j-j. per- 
mutes simultaneously all up and down spins, and 7 = ±1 
only since Pfa = i. 45 Together with P n Sf = -SfP n , it 
follows directly that 

{N c ;S*\§Z + §* Ne \N c ;S') 
= -(N c ;-S z \St + S z Nc \N c ;-S z )<xS z 

(as \N C ; S z ) is a spin doublet). Such matrix elements ap- 
pear in the z-component of the contraction in Eq. 4, and 
by spin isotropy an effective model of 2CK form results: 

# ff c = JC k,n c S ■ (§lo + sro) + S c l N S ■ (s L0 - s m ) (6) 

where S is a spin-i operator for the lowest chain dou- 
blet, defined by S z = £ s * \N C ; S Z )S Z (N C ; S z \ and S± = 
\N C ; ±^)(N C ; T^\- The effective couplings are given by 

Jk S n c = (N c ; +\\S Z + S z Nc \N c - + \)J K , (7a) 
Sf, Nc = (Nc, +\\Sl + S z Nc \N c - +\)8 K . (7b) 

Numerical evaluation of Eq. 7 for odd N c yields an AF 
effective coupling, J^ N > 0, which is renormalized with 
respect to the bare coupling and diminishes as the chain 
length increases, as shown in Table I. 



We consider first the case where the inter-impurity ex- 
change couplings J are sufficiently large that only the 
ground state of the isolated spin chain is relevant in con- 
structing the effective low-energy model upon coupling to 
the leads. As detailed in Sec. II B, this means in practice 
J > Tjf£ , with Tjf£ the scale for single-channel Kondo 
quenching of a terminal spin to lead a, arising in the "un- 
coupled" J = limit. The lowest state of an AF-coupled 
odd-membered spin-i chain is of course a spin doublet, 
the components of which we label as \N C ;S Z = ±5). All 
other states are at least 0(J/N c ) higher in energy. 

1. Effective 2CK model 

To leading order in Jkol , the low-energy model is then 
obtained simply by projecting into the reduced Hilbert 
space of the lowest doublet for a chain of length N c , using 

1n c =J2\ N ^ SZ )( N o' SZ \- ( 3 ) 
The resultant Hamiltonian H^f — iM c H^"l^ c follows as 

-^eff = |"foiiV c (Sl + Sjv c )ljV c • (SLO + Sijo) 

+5^1jv c (Si + SjvJ1jv c • (§L0 - s fl0 ) 

where J K = \{Jkl + Jkr) and S K = \(Jkl - Jkr), 
and we use the symmetry ljv c Siljv c = ljv c Sjv c liv c - 



TABLE I: Effective couplings 



N c 


J!i,N c 1 Jr 


1 


1 


3 


2 
3 


5 


0.51... 


7 


0.42... 



Hence, for sufficiently low temperatures T -C J/N c , 
the low-energy behavior of the system is equivalent to the 
single-spin 2CK model. 14 In the L/R mirror-symmetric 
case, Jkl = Jkr [5k = = 0), the stable T = 
fixed point (FP) is then the infrared 2CK FP. The lowest 
spin-i state of the impurity chain is thus overscreened 
by conduction electrons; giving rise to a residual entropy 
of Simp = |ln(2) (k B = 1), a hallmark of the NFL 2CK 
ground state. 14 Overscreening sets in below a character- 
istic scale T K CK , given from perturbative scaling 1 as 

T 2 K CK - J C K S N^M-1/PJ C K S N C ), (8) 

where p = 1/ (2D) is the lead density of states per orbital 
(assumed to be uniform) and 2D is the bandwidth. 

By contrast, when strict L/R mirror symmetry is bro- 
ken via distinct exchange couplings to the two leads, 
Jkl ^ Jkr (ie. 6 K ^ 0), the 2CK FP is destabilized 
in the full model Eq. 1. This behavior is of course 
well-known from the single-spin 2CK model with chan- 
nel anisotropy, 1-3,10 ' 14 where the impurity local moment 
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FIG. 2: Entropy 5 imp (T)/ ln(2) ds T/D for chains of length 
N c = 1, 3, 5, 7 (solid, dotted, dashed and dot-dashed lines re- 
spectively) in the large inter-impurity coupling regime. Shown 
for pj — 0.15 and pJn = 0.075 in the L/i?-symmetric limit 
8k — 0. Inset: scaling collapse onto the universal 2CK curve. 



is fully screened by conduction electrons in the more 
strongly coupled lead, and a Fermi liquid ground state 
results. For 8k > 0, under rcnormalization on reduc- 
tion of the temperature/energy scale, the system flows 
to strong coupling (SC) with the left lead (Jkl — > oo) 
while the right lead decouples {J KR O). 1 " 3 ' 10 ' 14 The 
stable T = FP thus depends on the sign of 6k, with 
'SC:L' describing the lowest energy behavior for 5k > 
while 'SC:R' is stable for 5k < 0. The mirror-symmetric 
case, 8k = 0, is as such the quantum critical point sepa- 
rating phases where a Kondo singlet forms in either the 
L or R lead. 

In the full model, effective single- channel Kondo 
screening characteristic of flow to the Fermi liquid FP 
in channel- asymmetric systems, sets in below a charac- 
teristic scale which can likewise be obtained from pertur- 
bative scaling: 1 

Tk° ~ JIn c *M-1/pJk,n c ), (9) 

where J KN = J'kn + \&1kn I 1S the effective coupling 
between the lowest chain doublet state and the more 
strongly-coupled lead. 



2. NRG results: symmetric case 

The above picture indicates that in the mirror symmet- 
ric case, 5k = 0, the lowest energy behavior for all odd 
chains should be that of the single-spin 2CK model, 14 but 
with renormalized effective couplings (Eq. 7 and Table I) 
and hence from Eq. 8 a reduced 2CK scale, T]f K . 

We now analyze the full model, Eq. 1, for odd N c = 
1, 3, 5, 7 using Wilson's NRG technique, 57-59 employing a 
complete basis set of the Wilson chain 60 to calculate the 
full density matrix. 60 ' 61 Calculations here are typically 
performed using an NRG discretization parameter A = 3, 
retaining the lowest N s — 4000 states per iteration. We 



FIG. 3: Spin susceptibility 4Txi mp (T) vs T/D for the same 
parameters as Fig. 2. Inset: scaling collapse of Xim P (T) itself 
onto the universal 2CK curve, with characteristic NFL form 
T K CK Xim AT) = aHT K OK /T). 

consider hrst the impurity contribution 58 ' 59 to thermody- 
namics (r2)i mp = (fi) — (&)o, with (Cl)o denoting a ther- 
mal average in the absence of the impurity chain. We 
focus in particular on the entropy, Si mp (T), and the uni- 
form spin susceptibility, Ximp 

(T) = ((5 2 ) 2 ) imp /T (here 
S z refers to the spin of the entire system) ; the tempera- 
ture dependences of which provide clear signatures of the 
underlying FPs reached under renormalization on pro- 
gressive reduction of the temperature/energy scale. 58 ' 59 

Fig. 2 shows representative NRG results for the T- 
dependence of the entropy, for odd chains with N c — 
1,3,5,7. At high temperatures the impurities are effec- 
tively uncoupled, so the chain contribution to the entropy 
is Si mp = AT c ln(2), as seen clearly from Fig. 2. On the 
scale of T <~ J/N c , all but the lowest chain doublet are 
projected out, the entropy then dropping as expected 
to ln(2) in all cases. Renormalization group (RG) flow 
to this local moment (LM) FP 1 ' 58 marks the regime of 
validity of the effective single-spin 2CK model, Eq. 6. 
The local moment is then overscreened 1,14 by symmet- 
ric coupling to the two leads on an exponentially-small 
scale, T K CK \ the T — residual entropy in all cases be- 
ing |ln(2). T]f K itself is determined in practice from 
Sim-p{T K CK ) = | ln(2) (suitably between the characteris- 
tic LM and 2CK FP values). The exponential reduction 
of the 2CK scale with increasing chain length (evident 
in Fig. 2) is expected from Eq. 8, which depends on N c 
through the effective coupling given in Table I. The inset 
shows the data rescaled in terms of T/T K CK . The low- 
temperature behavior of all odd chains collapse to the 
universal scaling curve for the single-spin 2CK model (ie. 
the N c = 1 case, solid line), likewise known from e.g. the 
Bethe ansatz solution, 2-6 confirming the mapping of the 
full model to Eq. 6. 

The above FP structure and energy scales also natu- 
rally show up in the magnetic susceptibility, as seen in 
Fig. 3. At the highest temperatures, Txi mp = jN c , as 
expected 19,58 for N c free spins. Flow to the LM FP on 
the scale T ~ J/N c is again clearly evident in the drop 
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FIG. 4: T = spectrum Dpk{<jj) vs u/D, for the same 
parameters as Fig. 2. Inset: vs ui/T]f K , showing collapse 
to the scaling spectrum for odd chains. Grey dotted line: 
\ow-\u\/T]f K <§C 1 asymptotic behavior Dpniy) = i[l — 

b(\uj\/T'k CK ) 1/2 ]; grey dashed line: high-|cj|/T| CK > 1 scal- 
ing behavior, Dp K {uj) = A/[\n 2 (\u}\/T^ CK ) + B}. 
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FIG. 5: Kondo temperature T^ CK /D vs 1/pJx for chains 
of length N c — 1,3,5,7 (circles, squares, diamonds and stars 
respectively, with lines as guide to the eye) and common ex- 
change coupling pj = 0.15. Inset: scaling collapse to common 
linear form when plotted vs l/pJ^Nc (given by Eq. 7). 



to Tx; mp = j, corresponding to the expected Curie law 
behavior. On the scale ~ T%P K , the susceptibility drops 



to T X u 



0, which remains its T = value. The inset 



shows TfP K Xim-p vs T/T^P K , showing scaling collapse 
to the universal single-spin 2CK curve, and demonstrat- 
ing the characteristic NFL logarithmic divergence 14 of 
the susceptibility Ximp itself on approaching the 2CK FP. 

We turn now to dynamics, in particular the low-energy 
Kondo resonance embodied in the spectrum Dpk{ljj) = 
DpK,a{u) — — TrpTlvn.[t a (iv)]; where t a (co) = t(cj) is the t- 
matrix 19 for the a = L/R lead (equivalent by symmetry 
for 8k = 0), and px = Np is the total lead density of 
states. Using equation of motion techniques 19 ' 66 yields 



■np T t a (uj) = -pJ Ka G a i(u) 
with i = 1 for a = L and i = N r for a — i?, where 



(10) 



G ai (uj) = ((Sr.f a0i + S?f a0f ; S+fU + ^ ot » w (H) 



J i J aOf' /uJ 



and {{A\B)) U is the Fourier transform of the retarded 

correlator «i(*i); B(i 2 )» = -ifl(ti-t 2 )<{A(ti),B(t 2 )}). 
The correlator in Eq. 11 can be calculated directly 
within NRG, 59-61 hence enabling access to the spectrum 
Dpx{oj). However, an alternative expression for the t- 
matrix can be obtained in the spirit of Rcf. 67, and in 
the wide flat-band case considered here is simply 



■np T t a {uj) = -i 



1 + 



G a0 (u>) 



(12) 



where G a0 (uj) = ({f a0(T ; fl a))u, is the Green function for 
the '0'-orbital of the a = L/R Wilson chain. 57 The quo- 
tient of correlators in Eq. 12 is found to improve greatly 
numerical accuracy, and is employed in the following. 



Fig. 4 shows the resultant spectrum Dpx(j-o) vs lu/D 
for chains of length iV c = 1,3,5,7 with the same param- 
eters as Figs. 2 and 3 (noting that pk{u) = Pk{— w) 
since the model, Eq. 1, is particle- hole symmetric). 
The low-energy form of each spectrum naturally re- 
flects RG flow in the vicinity of the 2CK FP, as studied 
also in a variety of different models which exhibit 2CK 
behavior. 13 ' 42 ' 45 ' 68-71 . In particular, for all odd chains, 
a half-unitarity value is seen to arise at the Fermi level, 
Dpx(uJ = 0) = |; and collapse to the universal single- 
spin 2CK scaling spectrum is clearly evident in the inset 
to Fig. 4. The leading low-frequency asymptotic behavior 
is D Pk {lu) = \[l - b(\uj\/T K CK ) 1 / 2 } (in marked contrast 
to the form [1 — a(|cj|/T^) 2 ] characteristic 19 of RG flow 
near a Fermi liquid FP); and with which the numerics 
agree well for lo <C T K CK . At high frequencies w » T K CK 
by contrast, the leading asymptotic behavior of the scal- 
ing spectrum is D Pk (oj) = A/[ln 2 (\oj\/T K CK ) + B], which 
behavior is common to other models in which spin-flip 
scattering processes are important at high energies, 72,73 
such as the single-channel Anderson or Kondo models. 

Finally, we consider the evolution of the 2CK scale it- 
self as the impurity-lead coupling is varied in the mirror- 
symmetric case, shown vs \/pJx in Fig. 5 for chains of 
length N c = 1, 3, 5, 7. An exponential dependence of the 
2CK scale on the impurity-lead coupling is expected from 
Eq. 8, and seen clearly in the main panel. The differing 
slopes reflect the renormalization of the bare Kondo cou- 
pling Jk —> J'kn W1 th increasing impurity chain length 
(Table I); collapse to common linear behavior being ob- 
served in the inset where the 2CK scales are plotted vs 
l/pJ[( N , establishing thereby quantitative agreement 
with Eqs. 7,8 and Table I, and hence the mapping to 
the effective 2CK model, Eq. 6. 
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3. NRG results: asymmetric case 

We turn now to the channel-asymmetric case, Jkl ^ 
Jkr (ie. 5k ^ 0). The effective model Eq. 6 should 
describe the low-energy behavior of all odd chains, so 
the rich physics of the asymmetric single-spin 2CK 
model 1-3 ' 10,14 is thus expected for T < J/N c . As dis- 
cussed above, breaking L/R mirror symmetry is a rel- 
evant perturbation 10 to the 2CK FP, so FL physics 
will arise generically on the lowest energy scales. In- 
deed, in the limit of maximal asymmetry Jk = 5k, 
the right lead is completely decoupled in Eq. 6; pris- 
tine single-channel Kondo (1CK) screening by the left 
lead then results below a single-channel scale TjF K . 
For 5k "C TjP K , however, RG flow in the vicinity 
of the 2CK FP strongly affects the behavior at higher 
temperatures/energies. 1-3,10 ' 14 This is indeed seen in 
Fig. 6(A), where the entropy Si mp (T) vs T/D is shown 
for a representative system with iV c = 3, pJ = 0.25, 
pJ K = 0.125, varying P 5 K = 10" 1 , 10~ 2 , 10~ 3 , 10~ 4 , 
10~ 5 , 10~ 6 , 10~ 7 [lines (a)-(g)], successively approach- 
ing the quantum critical point at the symmetric limit, 
5 K = 0: line (h). The scale T K CK (for which pJ K = 
p5k) is pT K CK w 10~ 3 here; and sets the scale for 
the crossover from 'large' to 'small' channel asymmetry 
(p5k > pT K CK and < pTjf K respectively). 

At the highest temperatures the three impurity spins 
are effectively free, yielding trivially a common entropy 
Simp = 31n(2). As T is lowered, all but the lowest trimer 
doublet state is projected out, heralding flow to the LM 
FP, with characteristic 1,58 entropy Si mp = ln(2). For 
large channel asymmetry [e.g. lines (a),(b)], RG flow is 
then directly to the SC:L FP: the impurity spin is fully 
screened 1-3 ' 10,14 by formation of a Kondo singlet with 
conduction electrons in the left lead (5k > 0) below 
an effective single-channel Kondo scale T^ c , and hence 
S iln p = for T < T K C . The Kondo scale itself is given by 
Eq. 9 in the large 5k regime, 74 with an effective Kondo 
coupling Jf, Nc =3 = 1{Jk + \5 k \) for N c = 3 (sec Table I). 

In the case of smaller channel asymmetry, p5k <C 
pTjf K , RG flow to the stable Fermi liquid SC:L FP oc- 
curs via the critical FP, which is of course the 2CK FP. 
The chain spin- ^ associated with the LM FP is then fully 
screened in a two- stage process (Fig. 6(A)). All such sys- 
tems flow first to the 2CK FP on a common scale T%F K , 
given by Eq. 8. The entropy thus drops to Si mp = | ln(2), 
symptomatic 14 of overscreening; before being quenched 
completely below a scale T <~ T* characterizing 1-3 ' 10,14 
the flow to the SC:L FP (with T* defined in practice by 
Sim P (T*) = i In 2). A clear S imp = \ ln(2) plateau is thus 
seen for lines (d)-(g) in Fig. 6(A), with T* diminishing 
rapidly as the transition is approached. The evolution of 
T* on varying 5k for different odd chains is itself studied 
in Fig. 8, below; and the result in the small- 1 5k | regime 
is a characteristic power-law decay, 
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FIG. 6: (A): Entropy 5 imp (T)/ ln(2) vs T/D on progressively 
approaching the transition, for N c = 3. Shown for fixed pJ = 
0.25 and p.J K = 0.125, varying pS K = 10 -1 , 10 -2 , 10 -3 , 10 -4 , 
10~ 5 , 10 -6 , 10 -7 [lines (a)-(g)], with the symmetric point 
P&k = shown as line (h). (B): Scaling collapse to the 1CK 
curve (circles) for strong channel asymmetry, pSx = 10 -1 
and 10~ 2 , with iV c = 1,3,5,7. (C): For small asymmetry 
(p5k = 10 -6 and 10~ 7 ), and N c = 1,3,5,7. Showing univer- 
sality in the LM— >2CK FP crossover, on rescaling in terms 
of T/T%f K ; compared directly to the symmetric 2CK scal- 
ing curve (circles). (D): Same data as (C), now rescaled in 
terms of T/T*; showing universal crossover from the 2CK FP 
(S im p/ In 2 = ±) to the stable SC:L FP (5 imp = 0). 



with exponent v = 2, and common amplitudes A on 
approaching the transition at 5k = from cither side 
(as guaranteed by symmetry). Eq. 13 generalizes the 
known result 2,3,14 for the channel-anisotropic single-spin 
2CK model, and is expected from the mapping of the full 
chain model (Eq. 1) onto the effective model, Eq. 6. 

We now turn to the scaling behavior of the entropy for 
chains of different length, as demonstrated by the three 
universal curves given in panels (B)-(D) of Fig. 6. First, 
in Fig. 6(B) for N c — 1,3,5,7, we show strongly asym- 
metric systems (p5k ^> pT K CK ) with p5k = 10 -1 and 
10~ 2 . The data clearly collapse to common scaling form 
when scaled in terms of T/T K C \ indicative of universal 
one-stage quenching from the LM to the SC:L FP. Re- 
sults for the single-channel spin- \ Kondo model are also 
shown (circles), confirming that the crossover is charac- 
terized by effective single- channel Kondo screening. 

The situation is more subtle for weakly asymmetric 
systems p5k <S pTjf K 1 where two-stage quenching oc- 
curs from the LM FP, through the 2CK FP, to the fully 
quenched SC:L FP. As now shown, each of these stages 
separately exhibit universal scaling, in terms of the two 
distinct low-energy scales T]f K and T* respectively. In 
Figs. 6(C), (D) for N c = 1,3,5,7, systems close to the 
transition are shown. To determine the full universal 
curves, it is of course essential to obtain good scale sep- 
aration of T]f K and T*: here T]f K /T* > 10 8 . 

In Fig. 6(C), results are rescaled in terms of T/T]f K . 
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FIG. 7: (A): Spectra DpK,h(u) (solid lines) and DpK,R,{oS) 
(dashed) vs oj/D on progressively approaching the transition 
for N c — 3, using the same parameters as in Fig. 6. (B): 
Collapse of Dpk,l{w) to the 1CK scaling spectrum (circles) 
in the strongly asymmetric limit, p8j< = 1CP 1 and 1CF 2 , for 
N c = 1,3,5,7. The dotted line shows the asymptotic low- 
\uj\/T k c FL behavior D Pk ,l{^) = 1 - a(\^\/T K c ) 2 . (C): For 
small asymmetry (pSx = 10~ 6 and 10 -7 for N c = 1,3,5,7), 
universal behavior arising on rescaling DpK, a {<->j) for both 
a — L/R in terms of oj/T k ck , as compared with the sym- 
metric 2CK scaling curve (circles). The dotted line shows the 
asymptotic (T* < |w| < T K CK ) NFL behavior Dp K , a {u) = 
|[1 - b(\u>\/TjP K ) 1/2 }. (D): Same data rescaled in terms 
of T/T*, showing the universal low-temperature crossover 
for Dpk,l{oj) (solid lines) and DpK,n{u) = 1 — Dpk,l{u) 
(dashed). The T* < |w| < T K CK asymptotic behavior 
Dpk,l{lo) — i[l + c(\oj\/T*)~ 1 '" 2 ] is shown as a dot-dashed 
line; while for |w| -C T* FL behavior results, DpK,h{u) = 
1 -d{\uj\/T*) 2 (dotted line). 

Collapse to the universal scaling curve for the symmetric 
single-spin 2CK model (shown separately, circles) is seen 
clearly; the crossover from the LM FP {S lrap = In 2) to 
the 2CK FP (S[ mp = ^ In 2) being as such determined 
by the 2CK scale T K CK . By contrast, the universality of 
the crossover from the unstable 2CK FP to the stable 
low-T SC:L FP with S imp = 0, is shown in Fig. 6(D). 
Here, on rescaling in terms of T/T* the data collapse to 
a universal form controlled by the low-energy scale T*, 
itself vanishing (Eq. 13) as the quantum critical point 
5k = is approached. 

The FP structure and energy scales naturally show up 
also in dynamical quantities, such as the scattering t- 
matrix, and hence the spectra DpK, a {^)- These are con- 
sidered in Fig. 7. In panel (A), again for N c = 3, we focus 
on Dp KtL {uj) (solid lines) and Dp K ^ R {uj) (dashed) for sys- 
tems with the same parameters as Fig. 6(A). Results for 
8k < are not shown, the a = L and R spectra simply 
being exchanged under the transformation 5k — 8k- 

N c — 3 chains with pSk = 10 _1 and 10~ 2 [strong 
channel-asymmetry, lines (a) and (b)], show a character- 
istic resonance in the left-channel spectrum, Dpk,l{w), 



FIG. 8: Fermi liquid crossover temperature Tfl/D vs pSk for 
chains of length N c — 1, 3, 5, 7 (circles, squares, diamonds and 
stars) and couplings pj = 0.25 and pJx = 0.125. For large 
channel asymmetry, Tfl = T K : good agreement between 
the data and Eq. 9 (dashed lines) is seen for each N c in the 
regime \pSk\ ^ 0.025. The inset shows the small asymmetry 
behavior on a log-log plot. Here Tfl = T*, and the quadratic 
behavior of Eq. 13 is shown as the solid lines, onto which data 
fall cleanly for p\5 K \ < 0.01. 



on the scale \ui\ ~ T^ c ; with the Fermi level value 
in particular being Dpk,l(0) = 1- This is the single- 
channel Kondo resonance; as is physically natural for 
these strongly channel-asymmetric cases because effec- 
tive single- channel Kondo screening is operative, 10 with 
the behavior thus expected to be that of the single- 
channel Kondo or Anderson models. 19 In the particle- 
hole symmetric limit of the latter, the Friedel sum 
rule 65 guarantees satisfaction of the unitarity limit 
(Dpk,l(0) = 1); and in the scaling regime \co\ <C J/N c 
one expects the entire one-channel scaling spectrum to 
be recovered. This is considered in panel (B), where 
results for N c = 1,3,5,7 and P 8 K = 1CT 1 and 1CT 2 
are shown, rescaled in terms of uj/T k c : essentially per- 
fect agreement is seen with the universal scaling spec- 
trum for the single-channel Kondo model (shown sepa- 
rately as circles). For T K C < \w\ < J/N c the char- 
acteristic Dp K ,h{u) ~ A/[\n 2 (\uj\/T K c ) + B] behavior 
typical of 'high' energy spin-flip scattering 72 ' 73 arises; 
while for |w| -C T K C canonical Fermi liquid behavior, 19 
Dpk,l{u) = 1 — a(\uj\/T^ c ) 2 , occurs as expected. 

Spectra for the right lead/channel, Dpk,r{oj) [dashed 
lines (a) and (b) of panel (A)], are similarly described by 
the leading <~ 1/ ln 2 (|cj|/T^ c ) logarithms at high ener- 
gies. However, the upward renormalization of the effec- 
tive Kondo coupling to the right lead — and hence RG 
flow towards the SC:R FP — is cut off at \uj\ - T| c , 
below which frequency the impurity chain local moment 
becomes fully screened by strong coupling to the left lead. 
Thus Dpk,r{oj) = for \oj\ < T K C , as observed directly 
from the NRG results in panel (A). 

We now turn to lines (e)-(g) of Fig. 7(A), for N c = 3 
systems with much smaller channel-asymmetry, pSk <C 
pT K CK . For both a = L and R, a clear half- unitary 
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plateau of Dp K>a (uj) ~ \ arises for T* < |w| < T|- CAr , 
indicative of RG flow near the 2CK FP. For \uj\ ~ T* , 
however, flow to the Fermi liquid SC:L FP occurs, such 
that Dpk,l{u = 0) = 1 and Dpk,r(co = 0) = are again 
satisfied. As was seen from the entropy (Fig. 6), there 
are two universal scales in this regime associated with 
the crossover from the LM FP to the 2CK FP [see panel 
(C)] and from the 2CK FP to the SC:L FP [panel (D)]. 

In panel (C) of Fig. 7, results are shown for systems of 
chain length N c — 1,3,5,7 and small channel asymme- 
try, p&K = 1CT 6 and 1CT 7 . Each is rescaled in terms of 
ui/T%P K 7 and collapse to the universal symmetric 2CK 
curve (circles) is seen in all cases — for both a = L 
and R spectra. In particular, for \w\ <C T^P K the 
characteristic NFL behavior is obtained, DpK. a {^) = 
\[1 - b{\u J \/T^ K ) 1 / 2 ] (dotted line). 

Finally, panel (D) shows the same data, but rescaled 
now in terms of \lo\/T*. Two universal spectra emerge: 
one for Dpk,l{u) and one for Dpk,r{u). The two scal- 
ing spectra are however found to be related simply by 
Dpk,r{u>) = 1 — Dpk,l(w), so we need consider only 
Dpk,l(u) (solid line). For ui ^> T* the relevant L/R 
symmetry-breaking operator dominates, 10 driving RG 
flow away from the 2CK FP. Since the scaling dimen- 
sion of this operator 10 is i, one expects Dpk,l(u) — 
i[f + cdwI/T*)- 1 / 2 ] (as indeed found, dot-dashed line). 
By contrast, for |u>| <C T*, irrelevant operators 10 ' 58 af- 
fect the RG flow in the vicinity of the stable FL FP, 
so one expects the leading low-|w|/T* asymptotics to be 
Dp K ±(J) = 1 - d(\uj\/T*) 2 (dotted line). Good agree- 
ment with the numerics is seen in both regimes and for 
both a = L/R spectra in Fig. 7(D). 

As the critical point is approached [\Sk\ — > 0; lines 
(a)— >(h) in Fig. 7], the spectra fold progressively onto the 
2CK spectrum itself [line (h)] down to lower and lower 
frequencies. The scale T* describing flow away from the 
2CK FP vanishes according to Eq. 13, as evident from 
the dynamics shown in Fig. 7 (or the thermodynamics in 
Fig. 6). In Fig. 8 the evolution of the low-energy scale as a 
function of channel asymmetry, pSx, is examined. Since 
Tft C is the lowest energy scale of the problem in the large- 
\6k \ regime, while T* is the lowest scale for small \Sk\, we 
consider the generic crossover scale Tfl in Fig. 8, defined 
in practice from the entropy via Si mp (Tp^) — |ln(2); 
which as such characterizes the flow to the ultimate sta- 
ble FL FP, and hence complete screening of the impu- 
rity spin. T-ph/D is shown vs pSx for chains of length 
N c = 1, 3, 5, 7. The dashed lines in the main panel show 
comparison to the perturbative result for TfP given in 
Eq. 9 (with a prefactor 0(1) for each N c adjusted to fit 
the numerics). The inset shows the same data on a log- 
log scale, demonstrating the quadratic decay of T*, given 
by Eq. 13 in the small-|^x| regime (solid lines). 



B. Weak inter-impurity coupling 

The perturbative derivation of the effective 2CK model 
in Sec. II A 1 is valid for sufficiently large inter-impurity 
exchange couplings, in which regard we note that any 
bare energy scale larger than the the exponentially-small 
universal scales T]f K or Tfl may be considered 'large'. 

However in the J = limit, the physical behavior 
is clearly very different. Here, each terminal spin- 1 ; 
impurity undergoes the standard single- channel Kondo 
effect 19 with its attached a = L/R lead below a temper- 
ature T ~ TjPJf, while the remaining impurities remain 
free down to T = 0. This scale is associated with the 
flow to the strong coupling (SC) FP, 19 and is given from 
perturbative scaling 19 by 

Tk C a ~ Dy/pJ Ka exp(-l/pJ Ka ). (14) 

The question addressed in this section is: what is the 
physics for odd chains with small but finite AF inter- 
impurity coupling, J < Tjf^l 

A physically intuitive picture for the simplest N c = 3 
case is depicted schematically in Fig. 9 (and discussed 
further in Sec. II B 1 below). Impurity T' forms a single- 
channel Kondo singlet with the left lead, and impurity 
'3' likewise with the right lead. A Fermi liquid descrip- 
tion then applies, and the remaining states of each lead 
act as an effective bath of non-interacting electrons that 
participate in the screening of impurity '2'. A residual 
AF exchange coupling, mediated via the Kondo singlets, 
once again yields an effective 2CK model. 

The above scenario is supported by the behavior of 
single-channel systems involving side-coupled quantum 
dots. 33-36 In the simplest example of a dot dimer, two- 
stage Kondo screening is operative in the small inter- 
dot coupling regime 34 : the dot connected to the lead 
undergoes a spin-i Kondo effect on the scale Tk,i, while 
residual AF coupling between the remaining dot and the 
lead gives rise to a second Kondo effect for T ~ T^ 2 (<C 
Tk,i), leading thereby to complete screening of both dots 
on the lowest energy scales. 3 Similar mechanisms have 
been advanced to describe the low-energy behavior of 
an asymmetric two-channel two-impurity Kondo model 24 
and a triple quantum dot ring structure. 43 

In the present context of odd impurity chains coupled 
to two leads, the most interesting behavior is expected in 
the i/i?-symmetric case. Here the effective coupling to 
left and right leads is also symmetric, and hence the 2CK 
FP must describe the low-energy behavior of the system. 



1. Effective 2CK model for N c = 3 

Before considering NRG calculations, we first derive 
the effective 2CK model for the simplest N c = 3 case; 
using perturbative techniques and scaling arguments, 
and exploiting the Wilson chain representation 57-59 (see 
Fig. 9) as natural within an RG framework. 
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FIG. 9: Schematic illustration for the N c = 3 chain with small 
inter-impurity coupling J < T]f^, as discussed in text. 



The Wilson chain for lead a = L/R is defined 57-59 
by dividing the band up into logarithmic intervals, 
{±DA~ n } (n = 0, 1, 2, ..), and then discretizing it by re- 
taining only the symmetric combination of states within 
each interval. This Hamiltonian is then tridiagonalized 
to obtain a linear chain form, with the impurity system 
coupled at one end. 57-59 The N c = 3 Hamiltonian may 
thus be written in a dimcnsionless form Hn = H + Hi, 



Hi = J(Si ■ S 2 + S 2 • S 3 ) 



N-l 



(15) 



+ (y-fanaf a(n+l)cr + J a(n+l)<T J ana) 



a, a n—0 



where s a rj is given in Eq. 2 and the Wilson chain op- 
erators P ana are obtained recursively using the Lanc- 
zos algorithm. 57 The rescaled dimensionless couplings are 
given by Jk<x = ^pJko/An and J = 2pJ/A^ (where 
A N = i(l + A-^A"^- 1 )/ 2 ). For a flat-band lead 
density of states, the tunnel-coupling between Wilson 
chain orbitals takes the form 57 t n = A( Ar - 1 )/ 2 A _ ™/ 2 ^„, 
with the = £„(A) - 0(1). The full Hamiltonian 
is then recovered in the N — > oo limit 57-59 via H = 
lim^oo {DA N H N }. 

We consider first the limit of strong impurity-lead cou- 
pling, J Ka ^> max(J',t„); so that H in Eq. 15 favours 
formation of a pair of singlet states between the terminal 
impurities and the 'O'-orbital of their attached lead. The 
ground state of H thus comprises a 1-'L0' singlet (we 
denote it by \s;L)) and a 3-'Rl' singlet (denoted |s;i?)), 
as shown schematically in Fig. 9. 

Hi now acts perturbatively, and we project onto the 
lowest state of H n using the unity operator for the re- 
duced Hilbert space, l s = \s; L)\s; R)(s; R\(s; L\. An ef- 
fective Hamiltonian may be obtained using the Brillouin- 
Wigner perturbation expansion, 75 Hn — E n + Hf s + 
Hff + Hf$+... Here, E = 1 S H 1 S = -\{Jkl + Jkr) is 
merely a constant shift in energy, while Hf s = l s Hil s = 
Hl follows as 

JV-l 

H L = X/ X] tn ijanaf a (n+l)a + fl(n+l)cr fanaj (16) 
cr,a n—1 

and corresponds to a pair of free Wilson chains, with the 
'O'-orbital of each removed. The second-order term, Hfj 



contributes only potential scattering 19 , here omitted for 
clarity. An effective coupling between impurity '2' and 
the 'LI' and 'Rl' orbitals is generated only to third-order 
in Hi, given 75 by Hffi = l s Hi{E - H )~ 1 PHi(E a - 
H )~ 1 PHil s with P = 1 — l s a projector. Combining 
this with Eq. 15, a rather lengthy calculation yields 



Hfll — pJlCL&2 ' SLl + PJi<R S 2 ■ Sri 



(17) 



(omitting RG irrelevant terms), with the effective cou- 
pling of impurity '2' to the a — L/R lead given by 



pJi<a — 



20t 2 J 



Ka 



> 0. 



(18) 



The effective model H^ = Hl + Hf{f, describes thereby 
the residual AF coupling between impurity '2' and the 
'LI' and 'Rl' orbitals of a pair of leads (with impurities 
'1' and '3' and lead orbitals 'L0' and 'R0' removed); see 
Fig. 9. This is a model of 2CK form. 

The above analysis presupposes the existence of the 
local singlet states \s;L) and \s;R). However, we note 
that for J < T}?k (as given by Eq. 14), RG flow is 
expected near a Fermi liquid-type FP, comprising single- 
channel strong coupling states in each lead, with a free, 
disconnected local moment on impurity '2' (and which 
'SC x SC x LM' FP is of course stable only at the point 
J = 0). Renormalization of the impurity-lead coupling 
JKa — > Jxa on successive reduction of the tempera- 
ture/energy scale naturally results in incipient formation 
of Kondo singlets between each terminal impurity and its 
attached lead below T ~ Tpfjf for the a = L/R channel, 
respectively. An effective 2CK model should then result 
via the mechanism described above, where the local sin- 
glet states \s;L) and \s;R) are now Kondo singlets. 

The central question then is: how does Eq. 18 flow 
under renormalization? Specifically, what is the effective 
coupling pjf a for T ~ Tjff ? 

To answer this, recall first that the effective 
temperature 57-59 within the RG framework is related 
to the iteration number/Wilson chain length, N, via 
T ~ A _JV / 2 . The operators for the Wilson chain orbitals 
also scale with N. In particular, operators for the 'O'- 
orbital of the Wilson chain scale as 57,58 f a o a <~ A~ N / 4 . 
Thus Jko. — > Jkol ~ Ar N / 2 , since the impurity-lead ex- 
change coupling is associated with a pair of 'O'-orbital 
fermionic operators. The key result is thus that the 
renormalizcd impurity-lead coupling Jk<x ~ Tjfjf for 
^ ~ — m accord with the physical expectation 

that disruption of the a = L/R Kondo singlet costs an 
energy OiTjf^). By contrast, the coupling between the 
impurities, J, is not associated with any chain opera- 
tors, and hence does not get renormalized with N. Fur- 
ther, as pointed out in Rcf. 76, once the 'O'-orbital of a 
Wilson chain has been frozen out (e.g. by formation of 
a Kondo singlet), the '1 '-orbital operators then scale as 
faXa ~ hr N ^. Thus the renormalized tunnel-coupling 
t Q ~ K- N ' A , so that for T ~ T\P£ , j%/J Ka remains 
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FIG. 10: (A): Si mp (T)/ln(2) vs T/D as the inter-impurity 
coupling is reduced, for the symmetric N c — 3 case (pSx = 0). 
Shown for fixed pJa = 0.15, varying pj = \{T]f K /D), where 
A = 10" n/1 ° and n = 0, 5, 7, 9, 11, 13 [lines (a)-(f)], with T, 
determined from a J = calculation. Behavior for T T/£ 
is described by Sf=°(T) = ln(2)+2S i 1 I ^f (T) (diamonds), with 
^irn^ (T) for a single-channel Kondo model with the same Jk- 
(B): Comparison of behavior for N c — 3,5,7 (solid, dotted 



1CK 
K 

1CK 



and dashed lines) with J = \T K CK . (C): Results for iV c 
3,5,7 (with n = 9 and 11): scaling collapse to the universal 
2CK curve (circles) is seen in all cases. 



0(1). From Eq. 18, the renormalized effective Kondo 
coupling at T ~ T K C £ can then be estimated to have the 

functional dependence pJ^ a ~ J I^W.a ■ 

For simplicity we focus now on the mirror symmetric 
case, where Jkc = Jk and Tjf^ = T K CK , whence one 
has the effective low-energy Hamiltonian for N c = 3 



H 



off 



ff L + /)4 tt Wc S.(s L1 + s ffl ), 
with S = S2 and the effective coupling 



pJ < K,N a 



J 



rplCK 

1 K 



x(N c ) 



(19) 



(20) 



valid for T < T K C . Determination of the constant 
x(N c ) is obviously beyond the scope of this analysis, al- 
though it can be deduced directly from NRG calculations 
as demonstrated in the next section. 

2CK physics is thus expected for T ~ T K CK (< 
T K CK ), as given from perturbative scaling 1 by 



T, 



2CK 



K 



rp\CK n reff 
1 K P J , 



K,N, 



exp(-l/ / 9j^ ff Ar J, 



(21) 



where the physical origin of the prcfactor Tjf K is simply 
that the effective bandwidth of the problem is already 
reduced to <~ Tjf K at the temperature T <** 
which the effective model, Eq. 19, is valid. 



Tjf K , below 



2. NRG results for odd chains 

The physical picture for the N c = 3 system is thus 
clear, and we now turn to NRG results for odd chains 



of length N c — 3, 5, 7 in the regime of weak coupling 
between the impurities. For accurate numerics, we found 
it necessary to retain N s — 4000, 6000 and 12000 states 
per iteration for N c — 3, 5 and 7, since higher-energy 
chain states remain important down to T ~ J < Tjf K . 

Fig. 10(A) shows S imp (T) vs T/D for the N c = 3 
case discussed explicitly above; for a common p.Jx and 
with pj = \{T}f K /D), where A = IO""/ 10 and n = 
0, 5, 7, 9, 11, 13 [lines (a)-(f)]. T}f K was itself determined 
from a J = calculation which, modulo a free spin on 
impurity 2, is equivalent to two separate single-channel 
Kondo models with the same Jk- At high T, the trivial 
S'imp = 31n(2) behavior expected for three free spins-i 
arises in all cases. Line (a) crosses directly to S^p = 
iln(2) on the scale J = T K CK w T K CK , characteriz- 
ing flow to the 2CK FP. By contrast, lines (b)-(f) flow 
first to the SC x SC x LM FP (5 imp = ln(2)). We also 
show for comparison S: 7 = p °(T) = ln(2) + 2S^(T) (dia- 
monds), where (T) is the entropy of a single-channel 
Kondo model 19 with the same Kondo coupling; S: J ^(T) 
thus describes the entire temperature-dependence of the 
entropy for J = 0. Lines (c)-(f) follow this curve per- 
fectly for T ^> T K CK , as expected from the single-channel 
Kondo screening of impurities '1' and '3'. An interme- 
diate Si m p = hi(2) plateau is thus observed, the single- 
channel T K CK remaining constant while the two-channel 
scale T K CK diminishes rapidly as J is decreased. RG 
flow thus persists in the vicinity of the SC x SC x LM 
FP over an extended T-range, but below T <~ T K CK all 
systems are described by the 2CK FP, with residual en- 
tropy S^p = |ln(2). 

We now comment on the generic behavior expected for 
impurity chains with A^ c > 3, which is a physically nat- 
ural extension of the N c = 3 case above. Following the 
'removal' of the terminal impurities through formation of 
single-channel Kondo sing lets for T ^T K CK , the remain- 
ing odd (N c — 2) impurities form a residual spin-i on the 
scale T ~ J (< T K CK ). This doublet state now feels an 
effective coupling to the two leads via the mechanism de- 
scribed in Sec. II B 1, but with a further renormalization 
of the effective Kondo exchange, as expected from the dis- 
cussion in Sec. II A in the regime of large inter-impurity 
coupling. Extension of the above analysis for N c = 3, 
which we do not give here, then leads us to expect (as 
tested below, Fig. 12) that the form Eq. 20 should hold 
for odd iV c > 3, with ratios x(N c + 2)/x(N c = 3) which 
are the same as those inferred from Table I but with two 
sites excluded from the spin-chain (reflecting quenching 
of the terminal spins '1' and 'iV c ' to form Kondo sin- 
glets); i.e. from Table I that x(N c = 5) w jx(N c = 3) 
and x(N c = 7) w 0.51x(A^ c = 3). 

2CK physics is thus is expected for all odd chains in 
the small inter- impurity coupling regime below T K CK , 
as given by Eqs. 21,20. This is confirmed in Fig. 10(B), 
where we consider chains of length A^ c = 3, 5, 7, tak- 
ing J = \T K CK as an illustrative example. Again, 
at the highest temperatures T 3> J, one obtains 
S'imp = iV c ln(2). For sufficiently large separation 
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FIG. 11: (A): Spectra Dpx(w) vs uj/D as the inter-impurity 
coupling is reduced, for the symmetric N c = 3 case (with 
same parameters as Fig. 10). Behavior at |w| S> T K K is 
described by Dp]f K (uj) (diamond points) for a single-channel 
Kondo model. (B): Comparison of N c = 3,5,7 (solid, dotted, 
dashed lines) with J = \T K CK . (C): Spectra vs uj/T k ck for 
N c = 3, 5, 7 with n = 9 and 11: universal scaling collapse 
is seen. For T K CK < |w| < T K CK the asymptotic behavior 
is Dp K (oj) = 1 - A/\\n 2 (\uj\/T K CK ) + B] (grey dashed line), 
while for |w| < T K CK NFL scaling is observed: Dp K {u>) = 
|[1 + b{\co\/T K CK ) 1/2 ] (grey dotted line). The entire curve is 
described by Dp K {u) = 1-Dp 2 £ K (lo) (circles), with pjf K (Lo) 
the scaling spectrum of the standard 2CK model. 



between T}f K and J, one expects the entropy to drop 



first to 5 imp = (N c - 2) ln(2) on the scale T^ CK (due 



to single-channel Kondo quenching 19 of the terminal 
impurities), followed by a further drop for T <~ J to the 
LM value S imp = ln(2); although in Fig. 10(B) Tjf K 
and J are comparable, so no distinct (N c — 2) ln(2) 
plateau arises. In all cases, however, S{ mp = \ ln(2) is 
seen below T^f K , characteristic of flow to the stable 
2CK FP; 14 with scales TjP K that evidently diminish 
with increasing N c , as expected qualitatively from the 
above discussion. Fig. 10(C) shows the scaling curve 
obtained when results are plotted vs T/Tfp K (with 
n = 9, 11 chosen to ensure good scale separation between 
T^f K and Tjf K ). The universal curve is precisely that 
of the standard 2CK model (circles). 

Dynamics are now considered briefly, Fig. 11 (A) show- 
ing the T = Dp K {uo) for iV c = 3 chains with the 
same parameters as in Fig. 10(A). All systems show 
RG flow in the vicinity of the Fermi liquid-type SC x 
SC x LM FP (reflecting single-channel Kondo screen- 
ing of the terminal impurities), and hence an incipient 
single-channel Kondo resonance 19 in each channel. For 
comparison, we show Dp]^ K (uj) for a standard single- 
channel Kondo model with the same Kondo coupling 
(diamonds), which recovers perfectly the spectral behav- 
ior for \oj\ ^> T 2 P K . In particular, for lines (d)-(f) 

2CK 



in the range Tj*- 



1CK 
K 



characteristic FL 



behavior 19 Dp K (u) = 1 - a{\oj\/Tjf K ) 2 arises, and thus 
the unitarity limit Dpk{u) = 1 is reached in this interme- 
diate energy window. For any finite J, T^P K is however 



FIG. 12: Evolution of the two-channel Kondo scale: T%P K /D 
vs pj for symmetric chains of length N c = 3, 5, 7 (squares, 
diamonds and stars), and Kondo coupling pJx = 0.15 (with 
5k = 0). For large inter-impurity coupling, T K CK is given by 
Eqs. 7,8 (indicated by arrows). For each chain, the maximum 
value T K CK « T K CK is obtained for J w 10T K CK . For J < 
T K CK , the Kondo scale diminishes rapidly; as seen also in 
the inset, where ln(T K CK ) vs T K CK /J is shown, confirming 
Eqs. 20 and 21. The slopes yield a;(iV c = 3) » 1/6, x(N c = 
5) w 1/9 and x(N c = 7) » 1/12. 



always finite, so ultimately RG flow to the stable 2CK 
FP yields D Pk {lo) = \ for |w| < Tjf K . 

Fig. 11(B) shows spectra for iV c = 3,5,7 chains with 
common J = \T]^ K , as in Fig. 10(B). These display the 
same qualitative behavior as for iV c = 3, with a single- 
channel Kondo resonance appearing at \co\ ~ T]p K , be- 



fore crossing over to the 2CK FP for \co\ ~ Tff K . Since 
T%P K diminishes with increasing chain length, while 
Tjf K remains fixed, apparent FL behavior consequently 
persists down to lower energies for the longer chains. 

For N c — 3,5,7 chains with the same couplings as in 
Fig. 10(C), Fig. 11(C) shows the spectrum arising when 
results are shown vs uj/T]f K : collapse to a single uni- 
versal scaling curve is seen clearly. One might naively 
expect to obtain the 2CK scaling spectrum here, since 
in this regime an effective 2CK model (Eq. 19) describes 
the system. However this model is only valid after the 
single-channel Kondo effect has already taken place in 
each lead, conferring a phase shift of ir/2 to the conduc- 
tion electrons. 19 In consequence, the scaling spectrum 
in the small J limit is Dpx{^) = 1 — Dp 2 ^ K (u), with 
Dp 2 ^ !K (ui) the scaling spectrum of the regular single-spin 
2CK model. Comparison with the latter (circles) con- 
firms this directly. 

Finally, in Fig. 12 we analyze the variation of the two- 
channel Kondo scale as a function of the inter-impurity 
coupling strength, pj, and the chain length, iV c . The 
T K CK / D in the large- J limit (indicated by arrows) are in 
accord with Eqs. 7,8. However, as pj is decreased, Tff K 
first increases (reaching its maximum of T]f K w T]f K 
for J w 10Tx K in each case), then diminishes very 
rapidly for J <C T]^ K . The behavior for small inter- 
impurity coupling is seen most clearly in the inset, where 
\n{T]f K /D) is shown vs T}f K /J. The linear behav- 
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FIG. 13: Schematic of the JV C = 3 trimer with impurities 
'1' and '3' modelled as Anderson-like correlated levels, each 
tunnel-coupled to both source (s) and drain (d) leads, with 
hybridizations r a7 (where a = L/R and 7 = s/d). The 
impurities are exchange-coupled to eliminate L -H- R charge 
transfer processes that destroy 2CK physics, but conductance 
through impurity 1(3) in channel a — L(R) can be measured. 



ior confirms Eqs. 20,21, with the slopes yielding x(N c = 
3) w 1/6, x(N c = 5) w 1/9 and x(7V c = 7) w 1/12, as 
consistent with the expectation discussed above. 



III. GATE VOLTAGE EFFECTS AND 
CONDUCTANCE 

Having analyzed in detail the Heisenberg chain model, 
Eq. 1, we now consider a variant in which the terminal 
impurities are treated as correlated levels (dots) , tunnel- 
coupled to their respective leads. Each lead a = L and R 
can be 'split' into source (s) and drain (d), allowing the 
conductance through dot '1' (or l N c ') to be measured; 
see Fig. 13. The Hamiltonian we study is "H Nc — Hl + 
Hc c + "H^ c b , where the four equivalent non-interacting 
leads are given by 



a— L/R k,(7 
f—s/d 



(22) 



and the impurity chain is described by 



JV C -1 

W^ c = J Sj • Sj+i 

i=l 



(23) 



+ U (ni T n u + h N ^n Nci ,) + e(hi + h Nc ), 



where hi = ^ a h ia = ^ CT d\ a d irT is the number operator 
for dot i = 1 or N c , e is its level energy and U its Coulomb 
repulsion/charging energy. In a real quantum dot device, 
the level energy is proportional to the gate voltage, e oc 
V g . The leads and chain are coupled via 



s/d 
k.cr 



(24) 

where V aj is the tunnel- coupling matrix element for the 
a = L/R and 7 = s/d lead. The hybridization strength 
follows as T Q7 = -kptV* (with px — Np the total lead 



density of states as before). Finally, a simple canonical 
transformation of the lead orbitals, via 



Cask<T — O'akcr COS^(7 Q 

) +5aka Sin(0 a ) 

Cadka =«ak ff sin(0 a ) - a akcr cos(6> Q ) 



(25) 



with tan(0 a ) = V a d/V as , yields an effective iwo-channel 
model with 



H L = J2 e k«lk ff a ak CT ( 26a ) 
a, k.cr 

<b = £ { V A k Ji, + V R a^d Nca + H.c.) (26b) 



where V 2 = V 2 S + V 2 d ; so that in particular for V 2 
Vl = V 2 (and hence T L = Y R = Y = np T V 2 ) the 
model is mirror-symmetric. To investigate 2CK physics 
on the lowest energy scales, this is the situation now 
considered. We also focus on the simplest example of 
the N c — 3 trimer, variants of which have been stud- 
ied recently in certain parameter regimes, 38-45 including 
exchange-coupled chain 42 and ring 43 ' 45 structures at half 
filling. As shown below, a physically intuitive pertur- 
bative treatment of the model for different fillings yield 
effective 2CK models - in spin and orbital sectors - from 
which the gate voltage-dependence of the 2CK scale can 
be identified. 



A. Effective low-energy models for N c — 3 

In the atomic limit (V — 0) of the isolated N c = 3 
trimer, the number of chain electrons jumps discontinu- 
ously between integer values N = n\ + n 3 + 1 = 1 — > 5 as 
the gate voltage V g oc e is varied (recall that impurity '2' 
is a strict spin-i). On tunnel-coupling to the leads, this 
Coulomb blockade (CB) staircase is naturally smoothed 
into a continuous crossover. Regimes of occupancy can 
still however be identified, and sufficiently deep within 
the CB valleys, A/" will be approximately integral. Here, 
a full <D{V 2 ) Schrieffer- Wolff (SW) transformation 19 ' 77 
can be performed in the strongly correlated regime of 
interest U » V, perturbatively eliminating virtual ex- 
citations into high-energy manifolds with (Af ± 1) chain 
electrons. 

In the atomic limit the ground state in any given TV- 
electron sector is a doublet, which we denote |A/";±i). 
Projecting onto the reduced (chain) Hilbert space of this 
doublet using the unity operator 



7=±t 



yields an effective model 



AT 



Hr 



-T/Cff 

rt A/\II> 



(27) 



where the 



leading 0(V 2 ) contribution arising from tunnel coupling 
to the leads (Eq. 26b with Vl = Vr and N c = 3, denoted 
as H') is given by the SW transformation 19 ' 77 



A/MI 



(28) 
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Here E Q = i^H^lM is the energy of the ground chain 
doublet (and retardation has as usual been neglected 19 ). 

In the following, we also exploit the particle-hole trans- 
formation 



5. 



~&2 ' 



a 



sz- 



,t 



(29) 



which yields directly TV — > (6— TV). The full Hamiltonian 
(parameterized by % = W{e) for given U, J, V) trans- 
forms as H{e) -> H(-e -U) + (2U - 4e). In general, the 
physical behavior of "H(e) and e — U) is equivalent, 
since the constant shift (2U— 4e) is irrelevant in the calcu- 
lation of observable quantities. Thus, the A/" = 1 f> 5 and 
TV = 2 o 4 sectors are related by the reflection about the 
particle-hole symmetric point e = —\U. Together with 
the singly-occupied TV — 3-elcctron case, three distinct 
regions of electron filling must in consequence arise. We 
now considered them in turn. 



2. TV = 2, 4 - electron regime: orbital 2CK effect 

As above, the TV = 2 and 4 electron regimes are related 
by the particle-hole transformation Eq. 29, so we consider 
explicitly only the TV = 4 case. The TV = 4 regime is 
the ground state of the free trimer over an e-interval of 
width J/2, specifically ~{U+\J) <e< -(U+\J). The 
ground state comprises a degenerate pair of spin singlets, 
since the spin-i on impurity '2' can form a local singlet 
with either '1' or '3' (the remaining site being doubly- 
occupied). Since the states are spin singlets, two-channel 
spm-Kondo physics will obviously not arise here. 

The TV = 4 states are however doubly degenerate, so 
may be associated with an orbital pseudospin (7~), and 
expressed as 



|TV = 4; V ) = % (« 2+CT ) t - SU\ 2+a)i 



dl 



) | vac), 



(33) 



1. TV = 3 - electron regime 

For —(U+jJ) <e< \ J, the atomic limit trimer 
ground state is a singly-occupied spin doublet 45 



|TV = 3; ^) = ^ [4 ff (4 r <-4/u) 

25. 



d\„dL 

1(7 OCT 



| vac), 



(30) 



where S z 



with a = ± for spins t/-|> an d £2 



4* is 



a spin raising/lowering operator. |vac) = J2a 2 \~ > CT2 > — ) 
defines the 'vacuum' state of the local (chain) Hilbert 
space, in which dots '1' and '3' are unoccupied, while '2' 
carries a free spin-i. 

Using Eq. 30 with Eqs. 27,28 leads eventually to the 
effective low energy model deep in the TV = 3 CB valley: 



UAf=3 — Hl + J'k^^S ■ (s L0 + sr ) 



(31) 



where we have omitted potential scattering contributions 
for clarity, and S is a spin- i operator for the lowest chain 

doublet, defined by S z = |TV = 3; S z )S z (Af = 3; S*\ 

and S ± = |TV = 3;±i)(TV = 3;t||- Eq. 31 is of 2CK 
form, with effective Kondo coupling 



+ 



1 



4e 
1 



5 J + 4e + 4U 5 J - 4e 



(32) 



which is AF throughout the entire TV w 3 sector. In 
the particle-hole symmetric Kondo limit in particular 
(e = — \U), one obtains pJ^ f A ^ =3 = \pJn to leading or- 
der in 1/U (with pJx = 8T/(nU) the effective exchange 
coupling of a single Anderson impurity 19 tunnel-coupled 
to leads); which as such is consistent with Eq. 7 and Ta- 
ble I for the N c = 3 Hcisenberg chain studied in Sec. II A. 

Two-channel Kondo physics thus arises in the TV w 3 
regime, with T^P K in particular given from Eqs. 8, 32. 



with T z = ^ for a = ±1. Projecting into this reduced 
Hilbert space using the TV = 4 unity operator Eq. 27 with 
the SW transformation Eq. 28, yields an effective orbital 
2CK model 



Hf =A =H L + Jf M=4 T ■ (f 0lT + r Qii ), (34) 
with an effective exchange coupling 



e + U-p 



+ 



e + U+p 



e + U+\j 



(35) 



which is AF within the TV w 4 regime. The trimer 
orbital pseudospin T~ is a spin-i operator defined by 

f z = T,T'W = ^T Z )T Z {N = A-T z \ and f± = 
|TV = 4;±i)(TV = 4; =F^ I . Similarly, we may define a 
lead pseudospin To jCT for each real spin a =^ / 4- as 

^0,<T — 2~(fl()crfL0cr ~ /iJOa /iJOa ) anQ ^ ^CKtr = /lOct/jJOct 

(with = (f + CT )t), where /2 0(T is given by Eq. 2b. 

The important 'pseudospin-flip' processes embodied in 
Eq. 34 correspond physically to moving an electron of 
given real spin from one lead to the other, while simulta- 
neously switching the trimer orbital participating in the 
local singlet, such that no net charge transfer occurs be- 
tween leads. Real spin a =t / I here plays the role of 
the channel index, and as such orbital 2CK physics is 
expected below T ~ Tjf K 1 as given by Eq. 8 using the 
effective coupling Eq. 35. 



3. TV = 1,5 - electron regime 

For e < —(U+jJ), the ground state in the atomic 
limit lies in the TV = 5-electron regime. It is a spin dou- 
blet, comprising the free spin-1/2 on 'impurity' 2, with 
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FIG. 14: Entropy S , imp (T) vs T/D for a trimer with tunnel- 
coupled leads. Shown for pj = 0.075, U/irT = 10, T/D = 
10~ 2 , varying e/irF — —5,-12.5,-14.5 (solid, dashed and 
dotted lines), respectively representing systems 'deep' in the 
Af = 3,4,5 CB valleys. Identical results are obtained for 
their particle-hole transformed counterparts. Inset: the low- 
temperature scaling behavior in terms of T/T% , compared 
with the standard 2CK model (circles). 

sites '1' and '3' each doubly occupied: \Af = 5;S Z ) = 
S% d[^d\^d^dl^\va,c) . Virtual excitations to the Af = 4- 
electron sectors are perturbatively eliminated to 0(V 2 ) 
by the SW transformation (Eq. 28), leading to 

Uf= 5 = H L + J C K \ =5 S 2 ■ (s L0 + §«,), (36) 
with effective coupling 

which is AF throughout the Af ~ 5-clcctron regime. 
Two-channel spin-Kondo physics is thus again expected 
below T ~ T$F K (given via Eq. 8). 

B. Thermodynamics and Scaling 

The physical picture is clear: 2CK physics, whether of 
spin- or orbital-character, is expected when sufficiently 
deep in each region of electron filling. We confirm this 
directly using NRG for the full trimer model in Fig. 14, 
where the entropy Si mp (T)/ ln(2) vs T/D is shown for 
fixed representative pJ, U/irT, and T/D, varying e/7rr 
for systems deep in the Af = 3, 4, 5-electron CB valleys. 

In each case, the high temperature T > U behavior is 
simply that of two free orbitals and a free spin, giving 
5 imp = 51n(2). The LM FP is reached directly as T is 
lowered, yielding S[ mp = ln(2); flow to the 2CK FP with 
characteristic 14 S imp = \ In(2) follows below T - T$F K . 
Upon rescaling in terms of T/TjP K (see inset), the sys- 
tems in each regime of filling collapse to the universal 



FIG. 15: Evolution of the two-channel Kondo scale Tjf K /D 
vs e/yrr for pj = 0.075, U/tvT = 10 and T/D = 10~ 2 . NRG 
results (points) are compared with Eq. 8 using effective Kondo 
couplings valid in the Af = 3 regime (Eq. 32, solid line), Af = 
2,4 regimes (Eq. 35, dashed lines) and the Af = 1, 5 regimes 
(Eq. 37, dotted lines). 

2CK curve (circles), thus confirming the effective low- 
energy models Eqs. 31, 34 and 36. 

Fig. 15 shows the evolution of the 2CK scale as the 
level energy (e oc V g ) is varied essentially continuously 
over a wide range of e/7rr, for systems with the same pJ, 
U/irT and T/D as in Fig. 14. NRG results (points) are 
compared with the perturbative result for T%f K given 
in Eq. 8, using the effective Kondo couplings valid in 
the M = 3 regime (Eq. 32), Af = 2,4 regimes (Eq. 35) 
and the Af = 1,5 regimes (Eq. 37). Throughout the 
majority of parameter space, the agreement is excellent; 
only at the boundary between regimes does the pertur- 
bative treatment (naturally) break down. Further, while 
the mechanism for overscreening changes from spin-2CK 
(odd-A/ - ) t° orbital-2CK (even-Af) across these bound- 
aries, T^P K itself is found from NRG to vary smoothly; 
the 2CK FP remaining the stable FP in all cases (includ- 
ing for e/7rr < -14 and > +4 in Fig. 15, where T$f K 
diminishes rapidly but nonetheless remains finite). 

C. Single- particle dynamics and Conductance 

We turn now to dynamics, focussing again on the 
spectrum — 7rpTlm[tz,(w)] = ttTDi(lo), where D\{ui) = 
— ilm[Gi(w)] with G\ the local retarded Green function 
for dot T'. We obtain it through the Dyson equation, 

[GiMrMGSMl^-EiM, (38) 

where G\{uj) is the non- interacting propagator (obtained 
for U = = J), and Ei(w) = Ef (w) - iEf(w) is the 
proper electron self-energy. The non- interacting G\ (ui) is 
simply 19 [GO(w)]- 1 =oj+-c-T(oj) (with lu+ = uj + i0+); 
where T(uj) = T^^-iT 1 (uj), with T 1 (to) = T (=nV 2 p T ) 
for all \u\ < D inside the band, and T r (lu = 0) = 0. 



An expression for Ei(a>) is readily obtained using equa- 
tion of motion methods, 19 ' 66 and is given by 

S 1 ( W ) = [G 1 H]- 1 {c/((rf 1T n u ;4 t )) w 

(39) 

+ij((rf u 4-+rf lt S|;dl t )) w J 

where the local Green function itself is Gi(u>) = 
{(d la ; d\ a )) u (independent of spin a in the absence of 
a magnetic field, and with G\(lo) — Gs(oj) in the mirror- 
symmetric systems considered). The self-energy can be 
calculated directly within the density matrix formulation 
of NRG 59 - 61 ' 67 viaEq. 39; with d(w) then obtained from 
Eq. 38. In particular, the local propagator for to = may 
be expressed simply as [Gi(0)] _1 = — e* +iT* in terms of 
the renormalized single-particle level e* and renormalized 
hybridization T* , given by 

e* = e + Ef (0) (40a) 
r*=r + E((0) (40b) 

in terms of the self-energy at cj = 0. The Fermi level 
value of the single-particle spectrum then follows as 

ttTD 1 (uj = 0) = — * (41) 

v ' r* i + (e*/r*) 2 v ' 

NRG results are considered in Fig. 16, panels (a)-(c) 
showing 7rr.Di(u;) vs uj/T for fillings Af = 3,4, 5, with the 
same parameters as used in Fig. 14. First, we comment 
briefly on the high-energy 'Hubbard satellites' clearly 
visible in D\{uj). As usual, 19 these reflect simple one- 
electron addition (co > 0) or subtraction (w < 0) from 
the isolated chain A/"-electron ground states. Their lo- 
cations, broadened somewhat on coupling to the leads, 
are thus readily understood from the V = atomic limit 
ground states of Sec. Ill A. Given their simplicity we 
do not comment further on them, save to note that in 
the Af = 3 particle-hole symmetric example of panel (a) , 
-Di(w) = Di(—u>) as expected; and that for the Af = 5 
example in panel (c), high-energy features are naturally 
observed only for uj < 0, corresponding to excitations to 
Af = 4-electron states. 

The most important feature of D\{uj) in Fig. 16 is of 
course the low-energy Kondo resonance, associated with 
RG flow in the vicinity of the stable 2CK FP. This is 
shown in panel (d) where spectra from the Af = 3, 4, 5 
regimes are again shown, but now rescaled in terms of 
uo/T]f K . Collapse to a single curve is seen, with the 
value at the Fermi level in particular pinned to ttTDi (lu — 
0) = 21 independent of e. The low-w asymptotics of the 
scaling spectrum (dashed line) are found to be 

M <<;1 

nTD^) T "~ ill-^M/TjF*) 172 ], (42) 

as consistent with behavior near the 2CK FP discussed 
in connection with a variety of related two-channel mod- 
els (see e.g. Refs. 13,42,45,68-71). Indeed, comparison 




uj/T oj/T£ ck 

FIG. 16: Left panels (a-c): T = single-particle spec- 
trum 7iTL>i(o;) vs oj/r. Shown for pj = 0.075, U/ttF = 
10, F/D = 10~ 2 , with e/7iT = -5 [Af = 3, panel (a)], 
-12.5 [Af « 4, panel (b)] and -14.5 [Af ~ 5, panel (c)j. 
Panel (d): collapse to the universal 2CK scaling spectrum 
(circles). The scaling behavior of the proper self-energy, 
plotted as E{(u;)/r and [Ef (cj) + e]/T, is shown in pan- 
els (e) and (f). Asymptotic |w| -C T]f K behavior de- 
scribed by ttPDiH = i[l - b{\ijj\/Tl CK ) 1/2 ] (dashed line); 
E(H/r = l+2&(|w|/T|° K ) 1/2 (dotted); and [Ef (w)+e]/r = 
-sgn(u;)26( \u |/lf cx ) 1/2 (dot-dashed) . 



to the spectrum Dpk(uj) for the 2CK model shows per- 
fect agreement in the low-energy scaling regime. This 
behavior is in striking contrast to that arising 19 in a FL 
phase for |w| < T K : ttTD(uj) = sin 2 (c5) - a x {u/T K ) - 
ci2(uj /Tk) 2 , describing the approach to the Fermi level 
value, which itself depends on the phase-shift, 8. For 
the Anderson model, 19 6 = ^n imp by the Friedel sum 
rule, 19 ' 65 with ni mp the 'excess' charge in the system in- 
duced by addition of the impurity. Thus, D(oj = 0) de- 
pends on the dot filling - and hence on the level energy e 
- in a regular FL. The situation is clearly quite different 
in the stable NFL phase obtained for the chain models 
studied in the present work; and we shall consider the 
analogue of the Friedel sum rule in Sec. HIE below. 

Further insight is however gained from the electron 
self-energy itself, the imaginary and real parts of which 
are shown respectively in panels (e) and (f) of Fig. 16. To 
emphasise the low-energy scaling of interest, wc show the 
results for Af = 3, 4, 5 in terms of u/T^ CK . The common 
asymptotic form for <C Tfp K is found to be 

S{(w)/r ~ l + 2b(\u;\/T 2 K CK ) 1 / 2 , (43a) 
EfH/r ~ -e/r-sgn^MM/Ti-^) 1 / 2 (43b) 

(with b ~ O(l) precisely the same constant as in Eq. 42). 
At the Fermi level in particular, E{(w = 0) = T, in con- 
trast to generic FL behavior E / (w = 0) = 0. Indeed, 
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extensive examination of NRG results over the entire pa- 
rameter space confirms Eq. 43 generally - for any value 
of the bare level energy e, and for all interaction strengths 
U/T and exchange couplings pj > 0. The renormalized 
level energy and hybridization then follow from Eq. 40 
as e* = and T* = 2T. From Eq. 41, the spectrum at 
the Fermi level is in consequence pinned to a universal 
half-unitarity value, itTDi(uj = 0) = \ for all underlying 
bare parameters, as illustrated in panel (d) of Fig. 16. 



1. Conductance 

To measure the differential conductance through dot 
'1', a bias voltage V s d is applied across the L source 
and drain leads, inducing a chemical potential difference 
A*s — — e Vsd- The L/R symmetry required to observe 
2CK physics on the lowest energy scales also requires of 
course that the same bias be applied across the R lead. 
Following Meir and Wingreen, 78 the zero-bias conduc- 
tance through dot '1' is given exactly by 



2e V 



dw(-)^M ttTD^lu) (44) 



where Di{uj) is the single-particle spectrum at equi- 
librium, /(w) = [e"/ T + I]" 1 , and r = T s + T d is 
the total hybridization as before. The dimcnsionless 
Go = 4T s rd/(T s + Td) 2 embodies simply the relative cou- 
pling to source and drain leads; such that for T s = Td 
Go = 1 is maximal, while in the extreme asymmetric 
limit T s >• Td (where the drain acts as a weak tunneling 
probe), G ~ 4r d /r s < 1. 

For T = 0, Eq. 44 reduces simply to 



2e 



G° = — G 7rrDi(w -0), 



(45) 



and Eq. 42 then gives a universal zero-bias conductance 



G°/G 



2 /h at T 



the gate voltage V g 



0, obtained for any value of 
This result is thus consis- 



tent with that known for related models in the singly- 
occupied Kondo limit which demonstrate 2CK behavior 
(see e.g. Refs. 15-18,24,42,45). For finite T, the Fermi 
level value of the spectrum has the same 
dependence 45 as the T = spectrum does of uj/T]f K 
(Eq. 42), viz. tiT£>i(w = 0;T) ~ ±[1 - b' {T/T 2 C CK ) 1 / 2 }. 
Combined with Eq. 42, Eq. 44 is then readily shown 45 to 
yield G C (T)/G ~ £ [1 - /T 2 K CK f'% with 7 = b' + 
by/nri( i ) and y/nr](^) ~ 1.07; which T-dependence is also 
known to arise for the single-spin 2CK model. 13,16 ' 17,70 

Calculating the conductance at finite bias is of course a 
different matter, and an exact (or even numerically exact) 
treatment of the underlying non-equilibrium physics is 
a formidable open problem. Here we merely make the 
simplifying approximation that the self-energy does not 
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eV sd /T 

FIG. 17: Conductance G c /(2G e 2 ft" 1 ) vs eV 3d /T for systems 
with same parameters as Figs. 14,16. Inset: scaling collapse 
to a universal curve with asymptotic eVsd <C T]f K behavior 
G c (V 3d ) = Goe 2 ^ 1 !! - b{eVsd/T 2 K CK ) 1/2 ] (dot-dashed line). 



depend explicitly on the bias voltage, 79 which leads to 

dfd^y 



G c — -tGo 
h 



ouj 



duj 



(46) 

where / 7 (w) = [e^ w ^-')/ T + l] 1 is the Fermi function for 
the 7 = s/d lead. While this approximation is exact both 
for V s d = 0, and for all V s d in the extreme asymmetric 
limit r s 3> T d 22 (conductance measured with a 'perfect 
STM'), in the standard case relevant to semiconductor 
quantum dot devices the leads are more symmetrically 
coupled. Here we consider a symmetric voltage split be- 
tween the leads, p s /d — ±^eV s d, for which Eq. 46 yields 



G c {V sd ) -J^GottT [Di(w - \eV sd ) + Di(w - -\eV ad )} 

for T = 0. This approximation thus allows us to work 
with single-particle spectra determined at equilibrium, 
obtained from a two-lead NRG calculation 59-61,67 as be- 
fore. 

Fig. 17 shows the resultant differential conductance 
G c /(2G e 2 h^ 1 ) vs eV sd /T, calculated using Eq. 47, for 
the same M = 3, 4, 5 systems as in Figs. 14 and 16. Since 
the conductance comprises a symmetrized combination of 
the total dot spectrum, similar features to that in Fig. 16 
are naturally observed. Peaks at high bias originate from 
the Hubbard satellites, and correspond to simple single- 
electron sequential tunelling processes. Importantly, the 
Kondo resonance also of course shows up, with the zero- 
bias value of G°/G = e 2 /h arising for eV sd < T]f K in 
each case. Universal scaling of the conductance in terms 
of eV s d/T^P K is also shown in the inset, demonstrating 
in particular the eV s d <C T]f K asymptotic form 



eV 3 , 



G c (V sd ) 



h 



Go[l-b{eV sd /T 2 K CK ) 1 / 2 ], (48) 



which behavior is likewise known 15,16 ' 18 in the NFL 
regime of the 2CK device constructed in Ref. 18. 
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D. Phase shifts and the S-matrix 

We now consider the leading low-w behavior of the 
single-particle scattering S-matrix, S(lj) — e 2lS ^\ and 
associated phase shift 8{uo) = S r (uj) + i8i(uj). The S- 
matrix is given by 65 

S(u) = l-2iYGi{tjj), (49a) 

with TG\(lo) related to the ti, (w)-matrix by 

irp T t L (Lj) = rGi(w) (49b) 

such that (as in Sec. IIIC) — 7rpTlmtL(o;) = ttTDi(uj). It 
follows from Eq. 49a that 



7rrDi(w) = | {l - e- 28 '^ cos(25*(w))} , 
leading in particular to the limiting behavior: 
sm 2 (S R (ujj) 



: S^uj) -> 
: Si(u) — > oo 



(50) 



(51) 



To obtain S(u>) it is convenient to express the propa- 
gator as G\(u) — [A(uj) + i-B(w)] -1 , where 

A(lu) = w-e-Ef(w), B(u) = T + S{(w) (52) 

such that A(0) = -e* and B(0) = V* by Eq. 40 (and for 
simplicity we take here the wide-band limit for the lead 
density of states, T(w) = — iT, which does not affect any 
of the following results). From Eq. 49a it follows that 



2iS(u) 



Aju) -iB'jbj) 
A(uj)+iB(uj) 



where 



B'(u) = r-E((w). 



(53) 



(54) 



First consider the familiar situation that would arise 
if the system was a regular Fermi liquid, for which 
E{(w = 0) = 0. In this case, B(0) = B'(0), and 
Eq. 53 yields 5(0) = arg[Gi(0)]. The S-matrix is then 
unitary at the Fermi level, ^(O)! 2 = 1; and, since 
Sj(0) = 0, the Fermi level spectrum follows from Eq. 51 
as tiT£>i(0) = sin 2 ((5 fl (0)). 

The situation is of course quite different for the present 
problem. The low-frequency behavior of the self-energy 
is given by Eqs. 43, and from which Eqs. 52-54 yield 



rp2CK 



[1 - % sgnH] (55) 



as the leading asymptotic form for \uj\/Tfp K — > 0; i.e. 

8r(uj) ~ iarctan[— sgn(w)] (56a) 



-25 1 (u) 



V2b 



T 2CKJ 



(56b) 
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FIG. 18: Luttinger integral Il/t vs level energy e/irF for sys- 
tems with pj = 0.075, U/ttF = 10 and T/D = 10~ 2 . Direct 
calculation via Eq. 63 shown as points, compared with Eq. 64 
(line) using ni mp determined from a standard thermodynamic 
NRG calculation. 



In evident contrast to a FL, the imaginary part of the 
phase shift thus diverges logarithmically as co — > 0, 



- i ln(M/T^), 



(57) 



the divergence itself reflecting (see Eq. 51) the pin- 
ning of the Fermi level spectrum to a half-unitary value 
(Sec. IIIC). In consequence, the S-matrix vanishes at 
the Fermi level, S(0) = 0, as known for the single spin-i 
2CK model. 80-82 This does not of course mean that an 
electron sent in to scatter off the dot is 'absorbed' (the 
conductance being generically non-zero), but rather that 
electrons scatter completely into collective excitations, 
characteristic of the NFL state. 80,82 

Notice also from Eq. 56a that the real part of the phase 
shift is discontinuous across the Fermi level, and that 
Sr(oj = 0±) 7^ arg[Gi(0)] (again in contrast to a FL). 
cos(2(5r(0)) = \ j\pl is however continuous across u = 
and, combined with Eq. 56b, Eq. 50 recovers precisely 
Eq. 42 for the low-w asymptotics of ttTDi(uj). 



E. Friedel-Luttinger sum rule 

We now consider further implications of the pinning 
of the T = Fermi level spectrum, ttTDi(0) = -|, re- 
gardless of bare model parameters and even when the 
dot occupancies change drastically on varying the bare 
level energy e. In particular, we obtain an analogue of 
the Friedel sum rule 19 ' 65 - a Friedel-Luttinger sum rule 62 
- relating the Fermi level spectrum to the 'excess' charge 
induced by addition of the impurity chain, 19 via the Lut- 
tinger integral. 63 ' 64 

To this end, consider first the excess charge n imp , de- 
fined as the difference in charge of the entire system with 



and without the trimeric impurity chain; 9 and also n[ 



imp' 
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defined correspondingly but with only the two terminal 
dots ('1' and '3') of the chain removed. Since impurity '2' 
is a strict spin, it follows trivially that n\ mp — n[ mp + 1. 
Using e.g. equation of motion methods, 19 ' 66 it is readily 
shown that 



4 

"imp = Im 

7T 



1 - 



ar(w) 



(58) 



(noting that sites '1' and '3' are equivalent by symme- 
try). In practice, as expected physically, n[ mp differs 
negligibly from the charge 2(hi) on the terminal dots, 
to which it reduces precisely in the wide flat-band limit 
where T(lu) = —iT is constant. 

Note next that Eq. 41 can be written as 



itTD 1 (lj = 0) = (r/r*)sin 2 (6») 



(59) 



with 9 = arctan(r*/e*) (= arg[Gi(0)]). Equivalently, 
using arg[Gi(o; = — oo)] = 0, 



f° d 
9 = lm dw— InGi(w). 

7-oo du 



(60) 



But from the definition of the propagator, Gi(u>) = — 
e — T(u>) — Ei(w)] -1 , it follows that 



Ai nGlM = _ GlM 



1 dT(w) 
du 



9Si(w) 
dui 



(61) 

The Friedel-Luttinger sum rule then follows directly from 
Eq. 60 as 



6 = 4™im P + II 
where the Luttinger integral 63 ' 64 

r ° .0Ei(w) 



r < 

II = Im / do; Gi(cj)- 

J — oo 



(62) 



(63) 



involves integration over all energy scales. 

Again consider briefly the situation that would arise if 
the system were a normal FL. In this case the Luttinger 
integral vanishes 63 ' 64 regardless of bare model parame- 
ters, and 9 = <5r(0) (as in Sec. HID). Eq. 62 then re- 
duces to a Friedel sum rule, 19 ' 65 relating the static phase 
shift to the excess charge. 

The present problem is not of course a Fermi liquid, 
and II does not vanish. The single-particle spectrum is 
however ubiquitously pinned at ttTDi(ui = 0) = ^ (i.e. 
T* = 2r and e* =0), whence 9 = n/2 regardless of bare 
parameters. 

Eq. 62 in this case thus becomes a sum rule relating 
the Luttinger integral to the excess charge: 

II = J(2-n| mp ) =J(3-n imp ) (64) 

Under the particle-hole transformation Eq. 29 it is easily 
shown that n[ mp -> 4 — n' imp (or equivalently 

"imp ' 



6 — «imp)- Hence II — > —II under the transformation, 
and in particular vanishes at the particle-hole symmetric 
point e = — U/2 where n- imp — 3 identically. 

The behavior of n; mp as e.g. the level energy e is var- 
ied, is naturally a smoothed/continuous version of the 
Coulomb blockade staircase arising in the atomic limit 
(where on varying e the total number of electrons in the 
free chain, N, jumps discontinuously between integer val- 
ues characteristic of each CB valley). Il itself will thus 
reflect that variation; and sufficiently deep in each CB 
valley, where n\ mp (~ A/") is close to integral, each regime 
may be loosely associated with its own value of II ■ 

The above discussion is exemplified clearly by Fig. 18, 
where the Luttinger integral Il/tt is shown vs the level 
energy e/nT for systems with common U/ttT = 10 and 
pj = 0.075. The points correspond to direct calculation 
of II via Eq. 63, using the full Green function and 
self-energy from NRG. The line is simply Eq. 64, using 
n; mp as determined from a standard thermodynamic 
NRG calculation. 58 The agreement is excellent over the 
wide range of electron fillings shown, the overall form 
of the curve reflecting the smoothed CB staircase as 
anticipated above. 

Our focus has been the N c = 3 trimer, but for longer 
(odd) chains one naturally expects the same 2CK physics 
to occur on low-energy scales. We have indeed confirmed 
this explicitly by NRG for the N c = 5 case. In particular, 
the single-particle spectrum of dot '1' at the Fermi level is 
again always pinned to half-unitarity, ttTDi(lj = 0) = | 
(whence the zero-bias conductance through a terminal 
dot remains G c /Gq = e 2 /h). The result 9 = n/2 thus 
holds generally, as does the Friedel-Luttinger sum rule 
Eq. 62, with n[ mp now related to the total excess charge 
by nimp = "!m P + (^c — 2); and in consequence the general 
result for the Luttinger integral for odd chains follows: 



j(2 - »!„p) 



lip) 



(65) 



IV. CONCLUDING REMARKS: REAL 
QUANTUM DOT SYSTEMS 

The exchange-coupled impurity chains studied in this 
paper may be considered as approximate low-energy 
models of quantum dot devices. In real systems, however, 
the dots are mutually tunnel-coupled rather than pure 
exchange-coupled; the 2CK fixed point is rendered un- 
stable by the inter-lead charge transfer that results, and 
the system crosses over to FL behavior on a low-energy 
scale T < Tfl- Experimental access to 2CK physics must 
thus contend with both channel anisotropy (as studied in 
Sec. II), and charge transfer terms. On the level of a toy 
model calculation, we now consider briefly the generic be- 
havior arising when the latter perturbation is included, 
considering explicitly the L/R mirror symmetric case (al- 
though the analysis is readily extended to include explicit 
channel anisotropy). 
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To motivate this, recall that a single (N c = 1) one- 
level quantum dot tunnel-coupled to two metallic leads 
does not of course exhibit 2CK physics. 20-22 This fol- 
lows from the Anderson Hamiltonian itself, i?And = 

Hdot + J2 a =L/R.l H L + #hyb]; wherC H dot = e(n dT + 

hdi) + Un dt h dl , HI = Ek,<r e k fl lk ff fl tt k ff and #hyb = 
<r( a ak(T^CT + H.c). Transforming canonically to 
even (e) and odd (o) lead orbitals 



ftektr =^75 (ctLkff + ClBka) 
a ok<T (flLkcr — Oflko-) • 



(66) 



i/ A nd is equivalent to H And = H dot + H e L + V2H^ yh , in 
which the dot couples solely to the e-lead; exhibiting as 
such single-channel physics only. 

In the singly-occupied dot regime, a low-energy spin- ^ 
Kondo model follows from a SW transformation 19,77 of 
HAnd, leading simply to -ffsw = H% + H^ + Hk with 



H K = JkS ■ (s LQ + s m ) + J LR S ■ s LR0 (67) 

(potential scattering is ignored); where s Q o is given by 
Eq. 2, and is defined as 



S.L.R0 - ^2 flou^crcr'faOcT' 



(68) 



with a = R, L for a = L, R. Eq. 67 consists formally of a 
symmetric 2CK model - the first term - together with a 
term S • s lro which transfers (cotunnels) charge between 
the leads. In fact, applying the transformation Eq. 66 
yields H sw = H L L + H° L + H K with 

H K = JkS ■ (s e0 + § o0 ) + JlrS ■ (§ e o - s o0 ), (69) 

where s e o and s D o are e/o lead spin densities; a model 
that is generically of channel-asymmetric 2CK form. But 
for the single-dot Anderson model itself the couplings 
are necessarily equal, Jlr = Jk- In this case the dot 
is exchange-coupled solely to the even lead, and Eq. 69 
reduces as it must to a single-channel Kondo model with 
Kondo coupling 2Jk- 

In systems comprising several tunnel-coupled quantum 
dots, however, cotunneling charge-transfer can be effec- 
tively suppressed, 24 with Jlr, <C Jk expected for longer 
chains (a simple estimate yielding Jlr/Jk ~ [t/U] N<: ~ 1 
with t the inter-dot tunnel coupling). Here we simply 
regard Eq. 67, with Jlr ^ Jk, as an effective toy model 
to mimic such effects in odd-iV c dot chains (with S rep- 
resenting the lowest chain doublet). From Eq. 69 it is 
clear that the resultant low-energy/T physics is then that 
of the channel-asymmetric 2CK model. 1 " 3 ' 10 ' 14 The 2CK 
FP is thus rendered unstable by the perturbation Jlr, 
any nascent 2CK state forming at T K CK being destroyed 
below the FL crossover scale Tfl; although inclusion of 
the Jlr term should not obscure the 2CK physics for 
T 3> Tfl, provided Jlr itself is sufficiently small. Note 
also that any direct inter-lead tunneling terms of the 
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FIG. 19: Upper panel: entropy 5 , i mp (T) vs T/D for a single- 
spin 2CK model with explicit left /right lead charge transfer 
(see Eq. 67). Shown for fixed pJx = 10 _1 , varying pJlr = 
KT 1 , 1(T 2 , 1CT 3 , 1(T 4 , 1(T 5 , 1(T 6 [lines (a)-(f)], successively 
approaching the pure 2CK limit with pJlr = 0, line (g). 
Lower panel: T = spectra Dpk(u) vs oj/D for the same 
systems. (T) and \Dp^ K (u) for a 1CK model with 

pJK = 0.2 (circles) are also shown for comparison. 



type Ek,k'(4 



Lk(T a Kk'cr 



H.c.) - as opposed to cotunnel- 
ing which intrinsically proceeds via the dot spin - are 
equivalent (through the transformation Eq. 66) to simple 
potential scattering in the even and odd channels. This 
does not destabilise the 2CK FP, 11 for which reason we 
do not include such terms here. 

The above scenario is explored in Fig. 19, where 
NRG results for Eq. 67 are shown. We fix pJk = 
1CT 1 , and vary the cotunneling term pJlr — 
1CT 1 , 1(T 2 , 1CT 3 , 1(T 4 , 1CT 5 and 1(T 6 [lines (a)-(f)], ap- 
proaching the pure 2CK limit Jlr. = [line (g)]. The 
top panel shows the entropy S , i mp (T)/ ln(2) vs T/D, 
from which the behavior associated with the channel- 
asymmetric 2CK model 1-3,10 ' 14 is seen to arise, as ex- 
pected from Eq. 69. In the extreme case Jk = Jlr. 
[line (a)], the odd channel is completely decoupled, no 
2CK physics occurs, and the behavior is that of a single- 
channel Kondo model (circles); the impurity entropy be- 
ing completely quenched below T ~ T K CK (cf. Eq. 14). 
For smaller J LR < T K CK however, the 2CK scale T K CK 
emerges (cf. Eq. 8), below which temperature a char- 
acteristic Si m p = |ln(2) plateau occurs [lines (d)-(g)]. 
Flow to the FL FP below T ~ Tfl, for all systems with 
Jlr 7^ 0, is then manifest in the final drop to Si mp = 0; 
the FL crossover scale Tfl being found to vanish as 



- FL 



A\J L r\ v :v = 2 (70) 
with exponent v = 2, just as expected from mapping to 
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the channel- asymmetric 2CK model, Eq. 69. 

We turn now to the T = spectra Dpx(u) = 
DpK,a{u) = — 7T/9Tlm[i Q (w)] for lead a = L or R (the 
two being equivalent by mirror symmetry) , shown in the 
lower panel of Fig. 19. The spectral behavior is rather 
different from what might naively be expected from the 
entropy, since no low-energy Tfl scale is apparent (com- 
pare e.g. to the channel-asymmetric 2CK models stud- 
ied in Fig. 7). This however reflects the fact that the 
model Eq. 67 is channel- asymmetric 2CK in the e/o basis 
(Eq. 69), rather than the L/R basis; and is readily under- 
stood using the transformation Eq. 66, from which one 
obtains DpK, a (u>) — \ [DpK, e {u) + DpK.o{^)] in terms of 
the spectra for e/o channels. For Jk = Jlr [line (a)], 
the odd channel is decoupled in Eq. 69, so DpK,u{u) = 
\Dp K ^ e {uj) = \Dp K JK (uj) 1 where Dp 1 K JK (uj) is the spec- 
trum for a 1CK model with Kondo coupling 2 Jk (as con- 
firmed explicitly (circles) in the lower panel of Fig. 19). 

However, for lines (c)-(f) (corresponding to Jlr <C 
T K CK ), the spectra are indistinguishable from the pure 
2CK spectrum, line (g) , over the entire range of frequen- 
cies. For T FL < \oo\ < T K CK , RG flow in the vicinity 
of the 2CK FP naturally results in the universal behav- 
for Dp K , e (u) = Dpk^uj) - i[l - b(\iu\/Tl CK y/\ as 
expected for a 2CK model with small even/odd channel 
asymmetry (see Fig. 7). Consequently Dpk(u) shows 
the same behavior. For |u>| <C Tfl by contrast, the spec- 



trum for the more strongly coupled e-channel has the 
asymptotic form DpK, e (<*>) = 1 — (^(M/Tfl) 2 while the 
weakly coupled o-channel is described by DpK,o{^>) = 
d(|u;|/TFL) 2 (sec Fig. 7 and discussion thereof). In- 
deed, we found in Sec. II A that the entire universal 
crossover to the Fermi liquid FP for the strongly cou- 
pled lead is related to that of the weakly coupled lead by 
DpK, e {u) = 1 — DpK,o(u). Thus, Dpx(u) — \ arises for 
all \u)\ -C T K CK , as indeed found. As such, the spectrum 
Dpk(uj) is effectively 'blind' to the Fermi liquid crossover 
induced by small finite Jlr- the 2CK FP appears to be 
stable on the lowest energy scales - although from e.g. 
the entropy we know that this is not the case. Ironi- 
cally, then, experiments that probe the t-matrix (such 
as measurement of the zero-bias conductance across dot 
'1') will always appear to yield 2CK physics, provided 
Jlr<T k ck . 
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